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a^ 

^jQ' In processes involving Coulomb-type initial- and final-state interactions, 

^ . the Gamow factor has been traditionally used to take into account these 

I additional interactions. The Gamow factor needs to be modified when the 

I magnitude of the effective coupling constant increases or when the veloc- 

' ity increases. For the production of a pair of particles under their mutual 
Coulomb-type interaction, we obtain the modification of the Gamow factor 

. in terms of the overlap of the Feynman amplitude with the relativistic wave 

i function of the two particles. As a first example, we study the modification 

I of the Gamow factor for the production of two bosons. The modification is 

' substantial when the coupling constant is large. 

a^ ■ 

0^ 



O ; I. INTRODUCTION 

Initial- and final- state interactions are important in many branches of tlieoretical physics 
involving the reaction or the production of particles. They have a great influence on the 
reaction rates or the production cross sections |]l|-p^. These initial- and final-state inter- 
', actions lead to a large enhancement of the cross section if the particles are subject to a 
strong attractive interaction; they can lead to a large suppression under a strong repulsive 
interaction. We shall use the term "the i^'-factor" to label the ratio of the cross section with 
the interaction to the corresponding quantity without the interaction. 

As is well known, for interactions such as the electric- Coulomb and color-Coulomb in- 
teraction V{r) = —a/r in non-relativistic physics, the effect of the initial- or final-state 
interactions leads to a A'-factor given by the Gamow- Sommerfeld factor The Gamow- 
Sommerfeld factor (or simply called, the Gamow factor) is given explicitly by 

= (1) 



X 



where 



V = -■ (2) 
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The magnitude of the relative velocity v is the ratio of the asymptotic momentum p to 
the energy in the relative coordinate system (see Eqs. ([T3|) and ( |T8|) below). Following 
Todorov [|1^], Crater and Van Alstine |T^, and Eqs. (21.13a)-(21.13c) of Crater et al. [|14| 



the relative velocity for the particles a and h is related to their center-of-mass energy ^/s by 



2711? 



(3) 



This gives v ~ 2yl — Am? / s when -^i ~ 2m and f — > 1 when s oo. This Gamow factor 
has been used to study initial- and final-state interactions in reaction processes. 

There are physical processes in which the coupling constant of the interaction between 
the particles can be quite large and the use of the Gamow factor to correct the initial- 
state and final-state interactions may not be adequate. For example, in the annihilation 
or the production of qq pairs, the interaction arising from the exchange of a gluon leads 
to a color-Coulomb interaction with a coupling constant a about 0.2-0.4, depending on the 
renormalization scale of the reaction process. Another example of strong coupling occurs 
in the case of a negatively charged particle in a nucleus with a large Z number. Such a 
large coupling constant will also lead to a modification of the Gamow factor as there are 
higher-order effects of the potential which are important when the coupling constant becomes 
large. One can mention, for example, the well-known case of the "Landau fall" which is the 
relativistic nonperturbative collapse of the wave function for an attractive Coulomb-type 



potential when the coupling constant exceeds a certain limit [0. Furthermore, although 
the effect of the interaction is very large for low relative velocities, it is useful to see how 
the effect varies as the velocity increases. 

While one sees the need to use the relativistic formalism to study the case with high 
relative velocities, one may wonder what special cases can be of interest to use a relativistic 
formalism for the case of low relative velocities. By the term "the relative velocity", we 
usually refer to the relative velocity between the particles in the asymptotic region of r — cx) 
where there is no interaction. However, when there is a strongly attractive interaction, the 
actual relative velocity depends on the spatial location. One can envisage that if the coupling 
constant is large, the motion of the two particles at small distances can become relativistic, 
even though the relative velocity at r — *■ cxo is small. Hence, it is necessary to use the 
relativistic formalism to study the effects of the mutual interaction with large coupling 
constants, even for the case of low asymptotic relative velocities at r oo. 

The i^'-factor for the Coulomb potential can be studied by examining the two-body 
wave function in the Klein-Gordon or the Dirac equation involving a Coulomb potential. 
Compared with the non-relativistic Schrodinger equation involving the Coulomb potential, 
there is an additional effective attractive potential, — |V^(r)p/2m^ (see Eq. (|16]) below), 
which leads to a non-trivial behavior when the coupling constant becomes large. In the case 
of fermions under the Coulomb interaction, there are further modifications associated with 
additional spin-dependent potential terms. 

The question of initial- and final-state interactions is also related to the question of the 
decay and the production of bound states when the interactions lead to the formation of 
bound states [^|T7|. Previously, the decay of the bound positronium ^5*0 state into two 



photons has been studied in the relativistic formalism by Crater |T^. In our present study. 



we are interested mainly in the case of two particles in the continuum. We wish to find out 
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how the mutual interaction may affect their reaction or production rates. 

The standard method of calculating the i^-factor is by evaluating the absolute square 
of the wave function at the origin. Such a method breaks down for the relativistic case as 



the wave function is infinite at the origin [|18[. The proper method to obtain the i^-factor 



is by taking the overlap of the relativistic wave function with the Feynman amplitude. As 
an illustration, we shall try out the method for the production of a pair of scalar particles 
interacting with a Coulomb-type final-state interaction. 



II. THE K-FACTOR 

The effects of the final- and initial-state interaction depend on the physical process. Here 
we shall be interested in the class of processes involving the reaction or the production of a 
pair of particles a and b, subject to the mutual interaction between a and b. For definiteness, 
we shall study the production process, as the A'-factor is the same for production or reaction. 

The simplest description of such a process is in terms of the perturbation theory which 
gives the amplitude for the production of this pair of particles a and b. The state of the 
ab pair after the reaction x + y—*a + b is represented by the state vector 

|$afe) = M{xy ^ a{P/2 + q)b{P/2 - q))\a{P/2 + q)b{P/2 - q)) (4) 

where A4{xy — > ab) is the Feynman amplitude for the x + y a + b process. For the 
two-particle system ab, we define the center-of-mass momentum P = a + b and the relative 
momentum q = {a — b)/2. 

On the other hand, under their mutual interaction which can be represented by a two- 
body potential V{r) between a and b, we can describe an ab pair with a center-of-mass 
momentum P as 

|vl>v^)=^(g)|P). (5) 

The probability amplitude for the production of an ab pair under their mutual interaction 
is obtained by taking the overlap of the amplitude in with the wave function in (j^). The 
overlap is the simplest in the ab center-of-mass system where P = {y/s, 0) and q = (0, q) 
and the (unnromalized) probability amplitude is ||16|,|l^,|19| 



(^v|$a.) = / ^^i^{ci)M{xy ^ a{q)b{-q)). (6) 



The A'-factor for the occurrence of ab in the state is then given by 

^_ \(^y\^ah)^ (Production cross section with final-state interaction) 
|(\l/o|^afe)P (Production cross section without final-state interaction) 

where |\l/o) is the state of the ab pair without their mutual interaction. We have used the 
unnormalized amplitude in Eq. (^ as any normalizsation constant will cancel out in the 
definition of the A'-factor in Eq. (|^). The result of the cross section calculated using the 
simple first-order diagram can be corrected to include the effects of the final-state interaction 
by multiplying with the A'-factor: 

/ Production cross section \ x ^ i P^^oduction cross section \ 
ywith final-state interaction j \^without final-state interaction j 
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III. KLEIN-GORDON EQUATION FOR THE COULOMB INTERACTION 



In this first study, in order to illustrate the main features of the effect and to avoid 
complications brought on by the spinor algebra, we shall carry out the procedures outlined 
above for the production of two scalar particles. 

We need to separate out the center-of-mass motion and the relative motion for these two 
particles. Consider first the case without a mutual interaction as in the region where the 
two particles are far apart. The two particles have 4- momenta pi and p2 and rest masses 



mi and m2- We introduce the total momentum P |ir3| , p 



P = Pl+P2 (9) 



and the relative momentum q 



where 



and 



q = wipi - W2P2 (10) 
s — ml + m? 

= — ir^ 

s — m? + 

W2 = 

2s 

with s = P^. We have the following identity 

pl-ml + pl-ml = {wj + wDp"^ + 2q^ -ml-ml = 0. (11) 

We can choose to work in the center-of-mass system in which P = {^/s, 0) and q = (0, q). 
The above equation can be written in terms of an effective energy e^, and a generalized 
reduced mass m^o as 

el-q^-ml = 0, (12) 

where 



s — mf — m2 



and 



(13) 



= — (14) 



Next, we study the system of masses mi = m2 = m interacting with a mutual Coulomb- 
type interaction 

V{r) = (15) 
r 

The equation of motion can be obtained from Eq. (|T2p by the canonical method of replacing 
e^; with e^, — V{r). The Klein-Gordon equation for the two-particle system under a mutual 
Coulomb-type interaction V{r) is 
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[e^-V{r)]'-q'-mljij{r) = 0. 
Writing ij^r) = Rni{r)Yim{0, 4>) in (p!6D, the equation for Rniij) is 

Rniir) = 0, 



d 2d in + l) 2ri o? 
\ H -A hi 

dz^ z dz z"^ z z"^ , 



where 2; = pr, p is the asymptotic momentum at r ^ 00 given by 



P 



(16) 



(17) 



(18) 



The wave function Rni{r) can be represented by the dimensionless variable z = pr and is 
characterized by two dimensionless parameters: rj = a/v and a^, where v = p/e^ is given 
by Eq. i). 

The solution of Eq. (|1^ is 



Rniir) 



Wia 



}}-e''''/\2izY''^/^e-"iFi{a, b, 2iz) 



where 



a = fi + - + ir], 



(19) 



(20) 



6 = 2/2 + 1, 



(21) 



(22) 



iFi is the confluent hypergeometric function, and the normalization constant has been de- 
termined by using the boundary condition that at r — > 00, Rni{r) sin(pr + 61) /pr with 
the Coulomb phase shift 61. For the S-state, the critical value of a is 1/2. 



IV. THE FEYNMAN AMPLITUDE AND THE OVERLAP WITH THE 

COULOMB WAVE FUNCTION 

To obtain the effect of the final-state interaction between a and b produced in a pair- 
production process, we consider the production of the pair of bosons from the fusion of 
two photons. Because the relevant factors associated with the mode of production will be 
cancelled out at the end in Eq. (|^), the results of the i^-factor depend only on the final-state 
interaction and is the same for a similar mode of production of the pair of bosons. 



ki{E,-k) 



HE, k) 



Pi{E, q) 
P2{E, -q) 



^ P2 



Fig. 1. Feynman diagrams included in the calculation. 
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The diagrams we include are shown in Fig. 1 which give the amphtude 

- iM{xy a{P/2 + q)h{P/2 - q)) = ie^ 



{2q + k) ■ El k ■ €2 ^ (2q - k) ■ €2 k ■ ei 



q — ky + w? 



(23) 



where k is the momentum of the photon, q is the momentum of one of the bosons, and 
ej is the polarization vector of the i-th photon. The overlap of the wave function with the 
Feynman amplitude is then 



{^v\^ab) = ie^ 



dq 



{2q + k) ■ El k ■ €2 {2q — k) ■ €2 k ■ ei 



{q + ky 



[q — ky + IV? 



(24) 



As the Coulomb wave function of Eq. ([T9|) is given in the configuration space, it is useful to 
write the above integral in terms of the wave function in configuration space. The latter is 
given by 



dq 



■ij{q)e 



iqr 



(25) 



In conventional applications, one expands the Feynman amplitude (^3]) in powers of q 
and keeps only the lowest order g-independent term J^q: 



M^Mo + 0{\q\). 



(26) 



In this approximation of taking only the leading term, Eqs. and ( P^ then give the usual 
i^-factor as the absolute square of the wave function iplr) at the origin 



K 



(r = 0)p. 



(27) 



However, such an approximation cannot be applied to our case with the relativistic wave 
function as the wave function, Eq. ([T9|), is infinite at the origin. To avoid this singular 
behavior, the full Feynman amplitude is needed to evaluate the overlap integral and the 
i^-factor in Eqs. (|2|) and (0). 

In terms of the spatial wave function, the overlap integral (|^) is 



—ikr 



{2ei ■ - - ei • fc}- 



47rr 



62 ■ k 



Akr 



{262- — + €2- fc}- 



A-nr 



ei ■ k 



(28) 



We shall specialize to the S-wave case with / = 0. Using the wave function of Eq. 
we carry out the above integration and obtain 



{^v I ^ab) = 2e\e^-k){e2-k] 



|r(a)| 

m 



n=0 



a)„r(| + /i + n) 



2ip 



n+fj.— ^ 



+ 



X 



2 3 /i + raS fi + n 5 2\ 2m^ + iJ, + n 5 fi + n 3 fi + n 5 2\ 

3 2 '4 2~''r^^^^^/WTW 2 '4 2~'2''^^ 



3 fi + n 1 _ fi + n 3 
4 ^ 2 '4 2 ' 2'^ ^ 



(29) 
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where 



6 = m + ip, 



(30) 



^ {m + ipy + k^' ^^^^ 

In Eq. ( p9|) (a)„ = a(a+ l)(a + 2)..(a + n — 1), with (a)o = 1. The quantity (6)„ is similarly 
defined. 



V. RESULTS FOR THE K-FACTOR 

We introduce the complex angle variable 

e = tan-i — ^ = T - 4 In (32) 

m + ip 4 4 K— p 

which is a relativistic measure of the relative motion between particles a and h. The real 
part of 9 is always 7r/4, and the imaginary part is negative, with a magnitude that is half of 
the rapidity of the produced particle in the center-of-mass system. 

In terms of the angle variable, Eq.(p9l) can be transformed as follows: 

{^v\^a,) = 2e^(ei • fc)(e2 ■ k)^^^e-^/'A, (33) 

where the factor A is 



g (a)„r(2z/)^ 2ip ^ 1 



X 



to ib)nn\ WFTP; A;2sin^ 

2 fsin(2z/-2)^ „ sm(2u - 1)9' 
' cos 9 



2z/-3 2u-2 2v 



2m 2v ^{mi(2v-\)9 ^ mi2v9\ . sm(2z/ - 1)^ 

+ sm 9 { — ^ — cos 9 ) - sm^ 9 — ^ — 

k 2v-2 \ 2v-\ 2v \ 2v-\ 



(34) 



and 1/ = 3/4 + (/i + n)/2. To obtain the i^-factor, we also need the overlap between the 
Feynman amplitude and the wave function without the final-state interaction. By using the 
wave function '?/'o('") = sin pr/pr for the S-state without the Coulomb potential, we obtained 
the amplitude for the case without the final-state interaction as given by 

(^o|$afe) = 2e2(ei-fc)(e2-fc)i3, (35) 

where the factor B is 

i3 = Jm{(cotr - 2m/A;)(r cotr - 1)} jpk, (36) 



and 9* is a complex conjugate of 9. Then the ratio between the absolute squares of Eqs.(^ 
and (|35|) is the relativistic expression of the i^-factor. 
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We can identify the factor |r(a)e'^''/^/r(6)p as closely related to the Gamow factor G{r]). 
One can show that 



T(a) 



m 



G{v) 



r(/i + i/2) 



r(2/i + 1) 



n 

j=Q \ ^ 



+J 



1 + 



:i+/x+2j)(i-/x) ' 



(38) 



Therefore, the proper treatment of the dynamics of the interacting particles leads to the 
modification of the Gamow factor G{ri) of Eq. (||) by a factor n given by 



K = G{7])k 



where 



K 



1 + 



(| + /x + 2j)(|-/x) ^ 



B 



(39) 



(40) 



and 



B 



E 

n=0 



2zp 



(a)«r(2z/) I 

(6), n! sin^ + A;^ 



■sin(2z/- 2)0 



2i/-3 



2t/-2 



cos 6* 



sin('2z/- r 



2m 2z/ , f sin(2z/ - 1)0 , sin2z/0' 
sm < cos - 



2z/-2""" [ 2z/- 1 2z/ 
/ |(A;/p)Xm{(A;cotr -2m)(rcotr - 1)}|^ 



2,sin 2z/ - 1 

sm'^^ ^ — 

2v-l 



2v-l 

2 



(41) 



In the limit of a; — or f ^ 0, the factor k goes to 1 and is consistent with the Gamow 
factor. 

We note that the center-of-mass energy yfs in units of the rest mass of the produced 
particle is a function of rj/a: 



m 



\ 



2 1 + 



7] /a 



(42) 



Various other kinematic variables, such as k/m = ^/s/2m and p/m = ^Js/4m'^ — 1 can be 
similarly expressed as a function of rija. From these relations and the relation between the 
i^'-factor and t] and a, we can study the i^-factor for the production of a pair of particles 
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in a specific kinematic configuration. 




12 3 4 

Fig. 2. The X-f actor versus 77 for various values of a. 

We sliow tlie beliavior of the i^-factor as a function of rj in Fig. 2 for various values of a. 
The solid curve gives the X-factor for a — 0.01 and the dotted curve gives the X-factor for 
a — 0.32. For a fixed value of a, the X-factor decreases as 77 decreases. This is consistent 
with the expectation that the effects of the final-state interaction diminish as the velocity 
becomes relativistic. The limiting value is K=l as rj = a/v — > a. Figure 2 also shows that 
for a given value of 77 = a/v, the /^-factor decreases as a increases. It should be noted that 
the same value of rj corresponds to different velocities v for different values of a. To see the 
effect of final-state interaction as a function of a for a fixed value of v, we plot in Fig. 3 
the X-factor as a function of v. As one observes, when the velocity is fixed, the i^-factor 
increases as the coupling constant increases, indicating a greater effect of the final-state 
interaction as a increases. For all values of a, the /^-factor decreases with v and goes to 
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unity as v approaches 1. The decrease is very rapid for small values of a. 
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Fig. 3. The /^-factor versus the velocity v for various 
values of a. 



It is of interest to see how the i^-factor obtained here is different from the Gamow factor 
in non-relativistic physics. In Fig. 4, we showed the ratio between the i^-factor and the 
Gamow factor for various values of a. As we expect, the ratio is almost 1 for weak coupling 
and the use of the Gamow factor is relatively safe there. However, if we increase a to 0.32, 
the ratio decreases significantly. The Gamow factor overestimates the magnitude of the final- 
state interaction. It cannot be used for the case with strong coupling. There is an effective 
screening of the long-range Coulomb interaction. As a consequence, the enhancement due 
to the long-range Coulomb-type interaction is reduced. It can also be observed in Fig. 4 
that the ratio of K/G{ri) is a relatively slowly varying function of r] for rj > 1 but drops 
down rapidly as rj decreases in the region of small rj. 

It is worth pointing out that the expansion of A in Eqs. (|3^) and (|41| ) is given as a 
series in powers of p/ \/5'^ + k"^ which increases as the velocity v increases. We still obtain 
convergent results for v up to about 0.8, but there is a limit on using such an expansion 
for greater velocities where p/V^^ + is too large to allow for a convergent term-by-term 
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summation. A different expansion method is needed. 
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Fig. 4. The ratio between the i^-factor and the Gamow 
factor G{ri) for various values of a. 



VI. CONCLUSIONS AND DISCUSSIONS 

The mutual final-state interaction between the produced particles has an effect on their 
rate of production. There will be similar effects if the particles interact via the initial-state 
interaction. The effects are simplest to be taken into account by using the method of the 
i^-factor. One calculates the rate for the process as though there were no initial- or final- 
state interactions, using, for example, the perturbation theory. The additional initial- or 
final-state interactions can be included by multiplying a ii'-factor as given by Eq. (^). 

For Coulomb-type interactions, the if-factor has been traditionally taken to be the 
Gamow factor obtained as the absolute square of the wave function at the origin of the 
relative coordinate. With relativistic Coulomb wave functions, the wave function at the 
origin is infinite and the usual method is not applicable. The /^-factor can be obtained as 
the overlap of the wave function with the Feynman amplitude. 

Our investigation of the i^'-factor for the case of the production of a pair of scalar particles 
indicates that there are substantial deviations from the Gamow factor when the strength 
of the coupling is large. In particular, the proper treatment reduces the magnitude of the 
Gamow factor significantly. The reason for this reduction is that in the pair production, 
there is an effective screening of the Coulomb-type interaction arising from the effective 
"exchange" of one of the produced particles. 
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We have presented an explicit formula for the relativistic modification of the Gamow 
factor for the production of a pair of bosons. Numerical results are also obtained to show 
the magnitude of the iiT-factor. The results of the i^-factor can be applied to a class of 
processes in which the boson particles are produced and interacting with a Coulomb-type 
interaction. 

The large modification of the Gamow factor for the production of two bosons studied 
here indicates the need to extend the present formalism to study the case of two fermions or 
two gluons. The application to fermions or gluons will be useful in the problem of production 
or reaction of quarks and gluons. As a pair of quarks or gluons interact with a strong color- 
Coulomb interaction with a coupling constant a about 0.2-0.4, the simple results from the 
present study indicate that the effects of the initial- or final-state interaction for quark and 
gluon will be large, and the i^-factor will be substantially different from what one obtains 
using the Gamow factor. We hope to study the modification of the Gamow factor for two 
produced fermions in our next investigation. 
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